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Abstract 

A parametrization of the lattice spacing (a) in terms of the bare coupling (/3) for the SU(3) 
Yang-Mills theory with the Wilson gauge action is given in a wide range of J. The Yang-Mills 
gradient flow with respect to the flow time t for the dimensionless observable, t-^P{E{t)), is 
utilized to determine the parametrization. With fine lattice spacings (6.3 < /3 < 7.5) and large 
lattice volumes (Ag = 64-128), the discretization and finite-volume errors are significantly reduced 
to the same level as the statistical error. 
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I. INTRODUCTION 


In the asymptotic-free gauge theories such as Yang-Mills theory and QCD, observables 
are obtained in the physical unit once the coupling constants of the theory are hxed at a 
chosen energy scale. In lattice gauge theory, this procedure is conveniently accomplished by 
identifying the energy scale with the inverse lattice spacing In actual applications, 

accurate determination of the relation between the lattice spacing and the gauge coupling is 
crucial not only for obtaining the observables in the physical unit but also for the continuum 
extrapolation of the numerical results. 

To derive such a relation, the reference scales, e.g. the string tension [1] and the Sommer 
scale [2] , have been widely used in the literature. In the determination of these quantities in 
the lattice unit, the heavy-quark potential V{r) for a certain range of the distance r needs 
to be htted, which inevitably introduces some systematic errors. 

Recently, a novel concept, the Yang-Mills gradient flow, was proposed in Ref. [3], and 
attracts wide attention both analytically and numerically. It provides us with conceptual 
and numerical advantages in lattice simulations [31116]. In particular, the gradient flow can 
be used to dehne proper energy-momentum tensor on the lattice [I7H22], which opens a 
new possibility to study the thermodynamics in lattice held theories [23l [2l]. To determine 
the reference scale with the gradient how, a dimensionless observable such as (see 

Eq. ([^ for the dehnition) is measured as a function of the how time t |3]. Then, t at which 
the observable takes a specihc value is used for the reference scale. This method does not 
require the htting procedure unlike the previous ones. Moreover, the statistical errors turned 
out to be substantially small compared to those of U(r) |3j. A variant of this method was 
also proposed in Ref. [5], where high-precision scale setting in lattice QCD is attempted. 

The purpose of the present paper is to determine the relation between the lattice spacing 
(a) and the bare coupling {/3 = Q/qq) in the SU(3) Yang-Mills theory with the Wilson 
gauge action over a wide range of (3 with high accuracy using the Yang-Mills gradient flow. 
Such a determination is useful for studying, e.g. the precision thermodynamics of the SU(3) 
Yang-Mills theory [23l [23]. In the previous studies with the Wilson gauge action [H l25] l26]. 
the range of [3 covered was 5.6 < (3 < 6.92. On the other hand, in the present study, by 
exploiting the beneht of the gradient flow together with the recent advance in computational 
power, we go into a weaker coupling region, f3 = 6.3-7.5. Moreover, to suppress the hnite 
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volume effect for the two finest lattices (/9 = 7.4 and 7.5), the lattice volume 128^ is taken 
for those cases. 

Following the procedures proposed in Refs. mi, we consider the dimensionless observ¬ 
ables and td{t^{E))/dt and £x their values to be X to determine the reference scale. 

We vary X in order to suppress both the lattice spacing and hnite volume effects below the 
magnitude of the statistical error. This enables us to derive accurate relation between a 
and (3 for the wide range of lattice spacing with an accuracy of less than 0.5%. 

This paper is organized as follows. In Sec. we introduce the gradient flow and define 
our reference scales. We then present our numerical results of the dimensionless observables 


in Sec. After describing our numerical setup, detailed analyses on the discretization and 
hnite volume effects are performed. The parametrization of the lattice spacing in terms of f3 
and comparisons with previous studies without the gradient how are also presented in this 


section. In Sec. IV, we give a brief summary. 


II. GRADIENT FLOW AND REFERENCE SCALES 


The gradient how is a continuous transformation of helds; for gauge helds, it is dehned 

by the diherential equation Ell], 

dA, 


dt 


_ „ 2 ^>S'YM(t) _ r, ^ 

i/o 


( 1 ) 


with the Yang-Mills action Symit) dehned by A^{t). Color indices are suppressed for sim¬ 
plicity. The A^(0) is identihed with the standard gauge held dehned in four dimensional 
space-time. The how time t, having a dimension of inverse mass-squared, controls the how 
into the extra dimension. The gauge held is transformed along the steepest descent direction 
of Syuit) as t increases. In the tree level, Eq. Q is rewritten as 


dAf^ 

dt 


= dydyA^ + (gauge dependent terms). 


( 2 ) 


which is essentially a dihusion equation. For positive t, therefore, the gradient how acts 
as the cooling of the gauge held with the smearing radius a/S^. In Ref. [6], it is rigorously 
proved that all composite operators composed of A^(t) take hnite valnes for t > 0. This 
property ensures that observables at t > 0 are automatically renormalized. 

One of the composite operators whose t dependence is extensively studied is 


m = 


(3) 


3 







where is the “field strength” composed of A^{t). The t dependence of the vacuum 

expectation value of Eq. ([^ for small t is obtained perturbatively up to next-to-leading order 

as 15 


e{E{t)) 


mi - mg? 

1287r2 


1 + ^9{qf + 0{g{qf) 


(4) 


with ki = A^"c( 117 e /3 -|- 52/9 — 31n3)/(47r) in the MS scheme and Nc being the number of 
color. The running coupling g{q) is defined at the scale of the smearing radius, q = I/a/s/- 
Equation Q implies that the dimensionless quantity is an increasing function of t for 

small flow time corresponding to the perturbative regime. As shown numerically in Ref. |1] 
on the lattice, t^{E{t)) is also a monotonically increasing function for large t corresponding 
to the non-perturbative regime. Therefore, the value of t at which t‘^{E{t)) takes a specific 
value X, i.e., the solution of the equation 


t^{E{t)) 


t=t 


X 


X, 


(5) 


is expected to be a unique dimensionful quantity, which can be used as a reference scale to 
introduce physical unit in lattice gauge theory. In Ref. |1], ^x=o .3 (sometimes called to) is 
used as the reference scale. In Ref. [5], a quantity defined by 




= .v 




( 6 ) 


is proposed as an alternative reference scale. In Ref. j5], a reference scale Wj^=o .3 (sometimes 
called Wo) is employed to set the scale and it was found that the discretization error of tco.s 
is suppressed more than that of to.s ia fall QCD [5]. 

In the present study we consider more general reference scales, and with X = 0.2, 
0.3 and 0.4. Larger X is preferable to suppress the lattice discretization error m, while the 
smearing radius would eventually hit the lattice boundary for too large X. We note that 
the numerical cost increases as X increases, since more time-steps are required for solving 
the differential equation Eq. ([^. We use wo .4 and wo .2 for the reference scales and introduce 
a new parametrization of the lattice spacing a in terms of the bare coupling /3 = 6/g^ hj a. 
hybrid use of these reference scales. 
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III. NUMERICAL RESULTS 


A. Simulation setup 


We perform numerical analyses of the SU(3) Yang-Mills theory with the Wilson plaquette 
action in the present study. Gauge conhgurations are generated by a combination of one 
heatbath and hve overrelaxation updates; we refer this set for updates as one Monte Carlo 
update. We perform 20,000 Monte Carlo updates for the thermalization from the cold start, 
and analyze conhgurations separated by 1, 000 updates after the thermalization. We take 
the periodic boundary condition with the lattice size N^. The values of (3 = 6/ Qq, Ng and 
the number of conhgurations Wonf are summarized in Table For f3 = 7.0, 7.2 and 7.4, we 
take two diherent values of Ng to investigate the hnite volume ehect. In our main analysis, 
we use the data with Ng = 64 for f3 = 6.3-6.9, Ng = 96 for jS = 7.0-7.2 and Ng = 128 
for f3 = 7.4-7.5. These data sets are indicated by * in the last column of the Table. As we 


will show in Sec. IIIC, the hnite volume ehect is well suppressed for these choices. 


Autocorrelation between diherent conhgurations is analyzed by the dependence of the 
jackknife statistical errors against the bin-size, Abm, together with the autocorrelation func¬ 
tion. For Ng = 64, we found no Abin dependence, so that diherent conhgurations are 
uncorrelated. For Ng = 96 and 128, the statistical error increases up to about 20% as 
Abin increases in Abin < 5. Then, we make the jackknife error estimate with Abin = 2 
in these cases. The autocorrelation function shows that the autocorrelation is not visible 
within statistics already with one separation of 1, 000 updates, which is consistent with the 
analysis of the bin-size dependence. 


The autocorrelation of the topological charge is known to become longer as the lattice 
spacing becomes hner due to the critical slowing down [3]. Our analyses with separations 
of 1, 000 updates may suher from this problem and therefore should be understood with 
reservation. New measurements of the topological charge with careful choice of the flow- 
time step size are left for our future work. 

The lattice discretization of the flow equation Eq. ([^ and that of the observable E are 
not unique [1]. We use the Wilson gauge action Syu for the flow equation in Eq. 0 H- 
To construct the operator A, we use the clover-type representation of unless otherwise 
stated. Other choices of Syu and E may reduce the discretization effect further m- To 
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solve the differential equation Eq. ([^ numerically, we use the second-order Runge-Kutta 
(RK) method. To estimate the numerical error of the RK method, we have solved the RK 
with two different integration step sizes, one of which is twice coarser than the other, for 
several configurations. By comparing these results we have checked that the numerical error 
of the RK method is within two orders of magnitude smaller than the statistical errors. 
(This procedure may be improved by the method suggested and tested in Refs. [8].) 

To illustrate the behaviors of t^{E{t)) and , we show these quantities with sta¬ 

tistical errors as a function of t/a? in Fig. [^for [5 = 6.3 and 6.7. To check the discretization 
effect on the operator, we compare t^{E{t)) and t-^t'^{E{t)) defined from the clover-type rep¬ 
resentation with those defined from E{t) = 2(1 — P{t)) using the average plaquette P{t) |1]. 
The figure shows that the difference between the two definitions is suppressed for large t. 
In particular, the difference is more suppressed in t^P{E{t)) than that in P{E(t)), i.e., the 
discretization effect in the former is smaller than the latter. The figure also shows that both 
quantities increase monotonically as t/a^ increases except for the small t/a? region where 
the lattice distortion effect is sizable. Furthermore, they show approximately linear behav¬ 
ior within the range 0.2 < X < 0.4. In Table the values of w^{a) and t^{a) obtained 
from Eqs. ([^ and (|^ with X = 0.2, 0.3 and 0.4 are given for each set of parameters. For 
simplicity, we will use the notations, and even for non-zero values of a. Note that the 
variances of and are approximately proportional to the inverse of physical volume of 
the lattice, while their dependence on a should be small because gauge invariant operators 
at positive flow time is automatically renormalized |1]. The statistical error shown in Table 
indeed shows such dependences on the spatial volume and lattice spacing. 


B. Estimation of for large X in high /3 region 

As pointed out in Ref. [5], the discretization error of Wq.s is smaller than to.s, so that we 
employ as the key reference scale in this paper. Note that the lattice artifact is expected 
to be smaller for larger X. 

In our simulations, we estimate wqa/o- for (3 = 7.4, 7.5 using the data at small flow time 
{ 100 , 2 /a) for these (3 and the extrapolation of the ratio wo, 4 ,/wo .2 to the continuum limit. We 
plot 9 data points for wo.a/wo .2 in the interval 6.3 < /3 < 7.2 as a function of in 

the left panel of Fig. In the figure, the results are shown for clover- and plaquette-type 
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TABLE I. Simulation parameters /3 and Ng, and number of configurations, Aconf- The result for 
the values of w^ja and y/t^ja are also presented. 


(3 

Ng 

-^conf 

woA/a 

wo.3/a 

WQ,2/a 

\/^0.4/o 

y/to^/a 

y/to^la 

Nga/wQ,2 


6.3 

64 

30 

3.208(7) 

2.877(5) 

2.460(4) 

3.269(4) 

2.835(3) 

2.254(2) 

26.02 


6.4 

64 

100 

3.697(5) 

3.317(4) 

2.837(3) 

3.765(3) 

3.263(3) 

2.590(2) 

22.56 

* 

6.5 

64 

49 

4.231(10) 

3.797(8) 

3.249(7) 

4.310(8) 

3.736(6) 

2.965(4) 

19.70 

* 

6.6 

64 

100 

4.857(11) 

4.356(9) 

3.725(7) 

4.938(8) 

4.277(6) 

3.390(4) 

17.18 

* 

6.7 

64 

30 

5.558(27) 

4.980(23) 

4.252(17) 

5.638(20) 

4.878(16) 

3.863(10) 

15.05 

* 

6.8 

64 

100 

6.300(20) 

5.652(17) 

4.833(13) 

6.406(16) 

5.548(12) 

4.395(08) 

13.24 

* 

6.9 

64 

30 

7.165(62) 

6.431(52) 

5.503(40) 

7.289(47) 

6.313(36) 

5.000(23) 

11.63 

* 

7.0 

64 

209 

8.177(34) 

7.322(28) 

6.250(21) 

8.287(25) 

7.168(19) 

5.674(12) 

10.24 


7.0 

96 

60 

8.137(21) 

7.297(18) 

6.236(13) 

8.264(16) 

7.154(12) 

5.665(08) 

15.39 

* 

7.2 

64 

204 

10.843(102) 

9.646(80) 

8.176(58) 

10.833(69) 

9.326(51) 

7.343(31) 

7.83 


7.2 

96 

53 

10.428(78) 

9.348(66) 

7.984(52) 

10.586(62) 

9.162(48) 

7.256(31) 

12.02 

* 

7.4 

96 

52 



10.426(106) 


11.927(98) 

9.415(64) 

9.21 


7.4 

128 

40 


12.084(61) 

10.306(42) 


11.808(40) 

9.337(28) 

12.42 

* 

7.5 128 

60 



11.706(72) 



10.601(42) 

10.94 

* 




FIG. 1. (the left panel) and (the right panel) as functions of t/a^ for (3 = 6.3 

and 6.7. Together with the result obtained with the clover-type operator for E, the one with E 
defined from the plaquette is also presented for each parameter. The statistical error is smaller 
than the width of the lines. 
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FIG. 2. Ratios wqa/wo .2 and \/tQA/'to .2 as a function of the lattice spacing. Squares represent 
the results with the clover-type representation of E, while circles show those with plaquette one. 
Solid and dashed lines represent 8-point and 9-point linear fits and the extrapolated values to the 
continuum limit are shown at a?‘jwQ 2 = 0 and —0.005, respectively. 


representations of E. The error bars of the data points are estimated by the jackknife method 
for the ratios and not for the individuals. To take the continuum limit, we perform a linear £t 
with all 9 points and that with 8 points by removing the data for the coarsest lattice. The 
same procedure for the continuum extrapolations will be adopted throughout this paper. 
The values in the continuum limit obtained from Fig. |^with the clover representation are 
(tCo. 4 /'W^o. 2 )a^.o = 1-3042(9) and 1.3047(10). The continuum extrapolation with the plaquette 
representation agrees well with this result. Similar extrapolation for shown in the right 
panel of Fig. [pleads to (\/to. 4 Ao. 2 )a-s>o = 1.4604(9) and 1.4600(11). 


The left panel of Fig. [^indicates that the lattice discretization error for the ratio wqa/wo .2 
is small. With this fact in mind, we estimate the values of Wqa/o. at /3 = 7.4 and 7.5 using 
the numerical results of tco. 2 /o and the linear fit shown in Fig. (left). The results are 
wqa/ol = 13.445(55)(5) at /d = 7.4 and wqa/ a = 15.272(95)(6) at = 7.5. The first error 
in the parenthesis is from the statistical error of wo, 2 /a and the fit parameters, while the 
second error is the systematic error obtained from the 9-point linear fit. The latter is more 
than one order of magnitude smaller than the former. 



















C. Effect of the finite voiume 


In order to investigate the finite volume effect, we have performed the numerical analyses 
for two different values of Ng for (3 = 7.0, 7.2 and 7.4 as shown in Table |T| where the spatial 
sizes in physical unit normalized by tco .25 -h/tco .2 = Nga/wo, 2 , are shown for each set of 
conhgurations. 

The comparison of the results with different Ng for f3 = 7.0 and 7.4 in Table shows that 
the values of w^/a and y/t^/a for different Ng agree within the statistics. On the other 
hand, the results for jS = 7.2 have statistically signihcant Ng dependence between Ng = 64 
and 96. These results suggest that the hnite volume effect modihes the numerical results 
for L/wo ,2 = 7.83, while the effect is not visible for L/wo ,2 > 9.21 in the present statistics. 
This is the reason why we use the data sets with * in the last column of Table |Tj 


D. Parametrization by the bare coupling f3 

For various practical applications, it is convenient to introduce a parametrization of the 
ratio wqa! 0 , in terms of jS. We have carried out such parametrization using four types of 
htting functions (polynomial type, one-loop type. Fade type, two-loop type) as summarized 
in Appendix All these htting functions can reproduce the numerical results in Table. 
well with three or four parameters. Among them, the three parameter £t motivated by 
the one-loop perturbation theory provides a reasonable result (y^/dof = 0.917) for 11 data 
points in 6.3 < (3 < 7.5 without over fitting: 

woi /47r2 37.422 143.84\ , , , . ^ 

-= exp I — 8.6853 H- - -— j [1 ± 0.004(stat.) ± 0.007(sys.)] (7) 

CL \ OO fJ jj J 

In the left panel of Fig. we show the numerical results of Table |T] normalized by the 
htting function Eq. ([^. The shaded band in Fig. |^is the error associated with the htting 
parameters in Eq. ([^. The results of some other htting functions in Appendix [A| normalized 
by Eq. ([^ are also plotted in Fig. They agree with each other within 0.5% in the range, 
6.3 < (3 < 7.5. In the right panel of Fig. the htting functions are plotted in the region 
beyond the present (3. Although the diherence among the curves grows as f3 becomes large, 
the deviation is still within 7% even aX (3 = 8.0. 
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{wq^I a) / (1-loop (3 param.)) 



1-loop (3 param.) 

1-loop (4 param.) 

polynomial (3 param.) - 

polynomial (4 param.) - 

Fade [1,1] — 

Fade [1,2] - 

Fade [2,1] 
lattice data '—B—' 


6.4 


6.6 


6.8 


(wq a)l (1-loop (3 param.)) 



0 


FIG. 3. (Left) Result of the three parameter fit of woa/cl given in Eq. 0. Squares are the data 
in Table |l] normalized by Eq. Q. Shaded band indicates the uncertainty originating from the errors 
of the fitted coefficients. Results with several other fitting functions normalized by Eq. Q are also 
plotted as well; see, Appendix [A} (Right) The same result is plotted for a wide range of 0. 




EIG. 4. Continuum extrapolations of the ratios y/t^/wQA and wqa/wqa- Squares represent the 
results with the clover-type representation of E, and circles are those with plaquette one. 

E. Continuum extrapolation of and 

We extract the continuum limits of \/t0/woA and w^/wqa with X = 0.2, 0.3 and 0.4 by 
plotting those quantities as a function of a?/ wq^ and making linear extrapolation to a = 0. 
Shown in Fig. |^are two examples of such extrapolation for X = 0.3. The resultant values 
are shown in Table [TT| where the statistical error in the hrst parenthesis is estimated by 
8-point linear extrapolation, while the systematic error in the second parenthesis is obtained 
by the difference between the 8-point and 9-point analyses as mentioned earlier. 
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y/t^/woA y/to^/wQA Wqa/wqa Wq.2/woa 
1.0164(32)(3) 0.8785(24)(0) 0.6952(18)(2) 0.8968(3)(2) 0.7665(6)(2) 

TABLE II. Continuum limits of y/t^/wQA and w^/wqa for X = 0.2, 0.3 and 0.4. The statistical 
and systematic errors are shown in the first and second parenthesis, respectively. 


rdwoA Tq/wqa d^WoA TcWqa tCo.4AMS 
1.328(21)(7) 2.587(45) 0.455(8) 0.285(5) 0.233(19) 

TABLE III. Relations of wqa with other reference scales in the continuum limit. For the details of 
the determination of these values and the estimate of the errors, see the text. 


F. Relation to other reference scales 


We now relate woa with other scales used in the literature. In Ref. [26], a scale Tc 
determined from the force F{r) between heavy quarks as dF{rc) = 0.65 is introduced and 
Tc/a was measured for four [3 values in the range 6.57 < (3 < 6.92. Using Eq. ([^, the 
result can be converted to the ratio rc/tco. 4 , which are plotted in Fig. The continuum 


extrapolation rc/tCo. 4 |a- 5 -o obtained from the hgure is given in the hrst column of Table III 


where the error in the second parenthesis includes the systematic error from the linear £t 
and also the uncertainly from the fitting function in Eq. ([^. 



FIG. 5. Continuum extrapolation of the ratio rdwQA using 3-point and 4-point linear hts. 
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The relation between wqa and the Sommer scale ro defined by r^FirQ) = 1.65 is also 
obtained by Tc/tq = 0.5133(24) given in Ref. |2S]; this is shown in the second column 


of Table III where the error takes into account all statistical and systematic ambiguities. 
The relation between tq and the string tension y/a is studied in Ref. [1] with the result 
= 1.178 within 1% uncertainty. The resulting value of s/awQA is shown in the third 


column of Table III We note that the continuum-extrapolated value of y/to/ro is estimated 
on coarser lattices in Refs. [2 m\. These values are consistent with our results within 
statistical errors. 


In Table III, we also show the relations of Wqa with the critical temperature of the 
deconfinement transition Tc and lambda parameter in the MS scheme, where we used 
= 0.625 ± 0.003(+0.004) [2H] and toAms = 0.602(48) [29]. We note that, in Ref. 


the value of /3 corresponding to the critical temperature with = 12 is obtained as /3c = 
6.3384. This together with Eq. ([^ leads to 


TcWqa = 0.2826(3), 


( 8 ) 


which is consistent with the value in Table IIIIl 

We note that Aj^ can be also determined by matching the tadpole improved coupling 
constant to that in the MS-scheme [30] . In Appendix we estimated the ratio tCo.4^MS 
by the same analysis as in Ref. [30] using our numerical data on plaquette and wqa/o,'. The 
result reads 

w;o.4Ams = 0.2388(5) (13), (9) 

which is consistent with the result in Table uni 


G. Relation to other parametrizations 

Let us compare our parametrization Eq. ([^ with those introduced in previous studies; 
Edwards, Heller and Klassen (5.6 < (3 < 6.5) [I], Alpha Collaboration (5.7 < (3 < 6.57) [20], 
Necco and Sommer (5.7 < (3 < 6.92) [20], and Diirr, Fodor, Hoelbling and Kurth (5.7 < 
f3 < 6.92) [31]. Each parametrization is based on the data obtained in the range of f3 given 
in the parentheses. In Refs. 123 and [20], fitting functions of the polynomial form are used, 
while in Refs. [1] and [3T] the fitting functions motivated by the perturbative formula are 
employed. 
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a) / (1-loop (3 param.)) 



FIG. 6. Comparison of the fitting function Eq. Q with previous results in Refs. [U [25l [26l [3T] . 
The numerical results on wqa/(i measured in the present study and Vc/a and ro/o in Ref. [26] are 
also plotted. Vc/a and r^/a are converted to rco. 4 -scale using the ratios in Table [nT| 


In Fig. we show the parametrizations of the above four references normalized by ours, 
Eq. ([^. For the conversion among different reference scales, we have used the ratios in Ta¬ 
ble |T^ The error of rc/wo. 4 , which dominates the ambiguity in the relations between wqa 
and other scales, is also shown by the shaded box in the hgure. The hgure shows that the 
parametrizations in the previous studies agree with ours within this error band in the range 
of /3, at which both htting functions are reliable. The parametrizations in Refs. [23] and [23] 
using polynomial ansatze have signihcant deviation from ours for /3 outside the applicable 
range, while those of the parametrizations in Refs. [1] and ina are much milder at the level 
of 4 — 6% for /3 = 7.5. 

In Fig. the numerical results on ro/a and Vc/a in Ref. [26] converted to wqa/ d using 


the values in Table |III| are also presented. The error bars of these points are the statistical 
errors in Ref. [26], and do not include the ones associated with rc/tco .4 and ro/tco. 4 - The 
hgure indicates that the r^/a and rda in Ref. [23] systematically deviate from our results 
as f3 becomes smaller. This may come from the discretization effect associated with the 
determination of tq on the lattice; these different parametrizations do not need to agree, 
because different determinations of the reference scales can differ by discretization effect. 
We also note that the statistical error in our numerical analysis of w^A/d is signihcantly 
smaller than the previous ones for 6.4 < (3 < 7.0. 
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IV. SUMMARY 


In this paper, we have performed an analysis of the flow time, t, dependence of 
and its logarithmic derivative for SU(3) Yang-Mills theory with 6.3 < (3 < 7.5 in large lattice 
volumes {Ng = 64-128). The results were utilized to parametrize the (3 dependence 
of the lattice spacing, Eq. Q. In our analysis, the reference scale is chosen to be tco. 4 , 
which is expected to suffer less discretization error than commonly used Wq.s and to. 3 - The 
discretization and hnite volume errors on our results are well suppressed with the present 
numerical settings. 

After the completion of this paper, we found the paper |32] in which Tc\/t^ = 0.2489(14) 
in SU(3) Yang-Mills theory is obtained. This is consistent with Tc\/t^ = 0.2483(7) obtained 
from the values in Table and Eq. (|^ . 
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Appendix A: Fit functions 

In this appendix, we show the four £t functions we employed to parametrize iCo. 4 /® in 
terms of (3'. 

1. Polynomials oi (3 — I3 q\ 

log { -) {13) = flo + Q‘i{I3 — /3o) + o.2{I3 — /3o)^ + n3(/3 — /3o)^- (^1) 

\ d / poly 
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fit func. 

xVdof 

polynomial (3 param.) 

1.496 

polynomial (4 param.) 

0.336 

1-loop (3 param.) 

0.917 

1-loop (4 param.) 

0.363 

[1,1] Fade (3 param.) 

0.377 

[2,1] Fade (4 param.) 

0.424 

[1,2] Fade (4 param.) 

0.424 


TABLE IV. Result of x^/dof of the fits for with various fitting functions Eqs. (Al)"(A3). 


2 . One-loop perturbation -|- polynomials of 1//?: 



47r^ 

lis" 


/^ + Co -1- 


/S 




+ ••• . 


(A 2 ) 


3. One-loop perturbation -|- Fade improved polynomials of l//d : 

1 4- di/j3 -\r <^ 2 //?^ 


log 




a /Pade^'^^ 331 4-ei//3- 1 -e 2 // 3 ^ 


(A3) 


4. Two-loop perturbation -|- polynomials of 1//3: 


log 


^ 0.4 \ , n 51 , ^ , C'i Co 

“ TAT + *^0 + ^ + ^ + 

d / 2-loop 00 Izl jj jj 


(A4) 


Here, Oj, Cj, di, e* and c' are fitting parameters. The polynomial form Eq. (Al) is the one 
used in Refs. [26] and [25] . In the three-parameter £t with Eqs. (Al), (A2) and (A4), we 


set Os, C 3 and Cg to zero, respectively. In the £t with Eq. (A3), we tried three parameter £t 
with d 2 = 62 = 0, four parameter fits with 62 = 0 and d 2 = 0. We refer each £t to as [1,1], 
[ 2 , 1 ] and [ 1 , 2 ], respectively. 

In Table IV we show y^/dof for each £t functions. For the polynomial £t, Eq. (Al), we 


take /3o = 7.0; the quality of the £t with other choices of (3o hardly changes. In the text, we 
employ the three parameter £t with Eq. (A2) to obtain Eq. ([^, since it provide reasonable 
X^/dof ~ 1 without over htting. 
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Appendix B: Determination of A-parameter 


In this appendix, we show the derivation of Eq.([^. In Ref. [30], ^o^ms analyzed with 
the data of ro/a in Ref. [20|. Here we adopt the same procedure by using the numerical 
results of wqa/o, and the average plaquette obtained in this study. Such an analysis allows 
us to determine the ratio tco. 4 Aiqs directly using the accurate data on hne lattices. 

The dimensionless parameter aAj^ can be obtained by matching the tadpole improved 
lattice perturbation theory. The boosted coupling constant gu is dehned by 

9h^9l{a)/ul (Bl) 

where Uq = P= (Tr t/n)/3. 

As for the choice of the renormalization scale and the running coupling constant, we take 
the following two methods: 


Method I 


at the scale 


and 


Method II 


with 




ap* = exp 


2bo 


£/□(«) \bo 


^ + (r^i -'>^ 2 ] aUa)+ 0{g^). 


a^MS = exp ( 1 > 


2br 


aAn = F^{gci{a)). 

This scheme corresponds to choosing a scale at 

t? \ F^ignia)) 


a/i= = exp 


2boJ F^^gaia) 


in Method I. 


(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 


In the 3-loop order, F^ {S = □, MS) is expressed as 
A ^ F%,(M)) = exp 


(B8) 
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where 


^ 6i 6? _ 

Pa = 


Ahl^hl-Ah^hl' 

In the [1,1] Fade approximation, it leads 


F[Ji](^/s(M)) = exp 


s hi h\- 2bobl 

PB = -jr2 + 


Abl^bj-Abobf 


26o^s 


bogl 




^1 

•2bo 


In SU(3) Yang-Mills theory, the coefficients are given by 


bn = 


11 

(47r)2 ’ 


bi = 


- M 2 > (Bll) 


with 


102 ^ _ 1 2857 

(47r)*’ ^ ^ (4ir)8 2 ’ 

t? = 0.1348680, 4° = 0.0217565. 


(B9) 


(BIO) 


(B12) 


We adopt Method II with Fade improvement to estimate the central value of t(;o. 4 Aj^, and 
use the results of Methods I and II without Fade improvement to estimate the systematic 
error. The results are summarized in Table IVl The values in the continuum limit are obtained 
by a linear £t as a function of a^/tCo ,4 without the coarsest lattice data (see Fig. [^. Then 
we hnd 

= 0.2388(5) (13) (B13) 

with the statistical and systematic errors. 
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(3 Ns plaquette 

woA/a 

W^0.4A]y[g 

Method I 

Method II Method II Bade 

6.3 64 0.622 420 85(30) 

3.208( 7) 0.2248(4) 

0.2253(4) 

0.2234(4) 

6.4 64 0.630 632 88(13) 

3.697( 5) 0.2280(3) 

0.2285(3) 

0.2266(3) 

6.5 64 0.638 361 33(35) 

4.231(10) 0.2298(5) 

0.2302(5) 

0.2285(5) 

6.6 64 0.645 669 58(12) 

4.857(11) 0.2327(5) 

0.2331(5) 

0.2314(5) 

6.7 64 0.652 608 39(39) 

5.558(27) 0.2351(11) 

0.2354(11) 

0.2338(11) 

6.8 64 0.659 215 11(11) 

6.300(20) 0.2354(7) 

0.2358(8) 

0.2342(7) 

6.9 64 0.665 522 54(33) 

7.165(62) 0.2367(20) 

0.2370(20) 

0.2355(20) 

7.0 96 0.671 556 729(89) 

8.137(21) 0.2378(6) 

0.2381(6) 

0.2367(6) 

7.2 96 0.682 891 86(22) 

10.428(78) 0.2389(18) 

0.2392(18) 

0.2379(18) 

oo 1 

oo 0.2399(5) 

0.2401(5) 

0.2388(5) 


TABLE V. Simulation parameters j3 and Ag- The plaquette value, wqa/cl, and rt;o. 4 Aj^ using 
Method I, II and II with Bade approximation. The last row corresponds to the values at the 
continuum limit obtained from linear extrapolation without using the coarsest lattice data at /3 = 
6.3 (the italic number). 



FIG. 7. Values of mo. 4 A^ estimated by the Method II with Bade improvement. The continuum 
limit is shown at a? ^ = 0. 
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